A037  754  AIR  FORCE  WEAPONS  LAB  KIRTLAND  AFB  N MEX  F/6  20/13 

THERMAL  RADIATIVE  LOSSES  FROM  GRAPHITE  CALORIMETERS  USED  IN  PUL— ETC (U) 
NOV  76  J A MACFARLANE 

UNCLASSIFIED  AFWL-TR-76-180  NL 


I * I 

1%37754 

El 

1 

L_ 

- 

END 

DATE 

FILMED 

4-77 

J 


n r 


AFWL-TR-76-1 80 


' ; 

fi 


THERMAL  RADIATIVE  LOSSES  FROM  GRAPHITE 
CALORIMETERS  USED  IN  PULSED  ELECTRON  BEAM 
EXPERIMENTS 


November  1976 


Final  Report 


€ ** 


\\  ■ ^ 


Approved  for  public  release;  distribution  unlimited. 


This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under 
Subtask  N99QAXAA121,  Work  Unit  10,  "Diagnostics  Development." 


4 ,f  r*  Um 

pi»PY  RVfc 

,A,^vy‘,wu 


Prepared  for 
Oirector 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 


w- 1 


JL  ** 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


AFWL-TR-  76- 1 80* 


This  final  report  was  prepared  by  the  Air  Force  Weapons  Laboratory, 

Kirtland  Air  Force  Base,  New  Mexico,  under  Job  Order  WDNA4501 . Lt  John 
A.  MacFarlane  ( DYV ) was  the  Laboratory  Project  Officer-in-Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a definitely  related  Government  procurement  operation, 
the  Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever, 
and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 


Chief,  Materials  Vulnerability/  Chief,  Technology  Division 

Hardening  Branch 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flFhon  Dal  a Entered) 


T"  AFWL-TR- 76-180 


REPORT  DOCUMENTATION  PAGE 

IfRgPOWT-mTMaPB 


2.  GOVT  ACCESSION  NO 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


' 


TITL&.14PH  Subtitle 1 


i i THERMAL  RADIATIVE  10SSES  FROM  GRAPHITE  CALORIM-  ( 

* Iters  used  in  j>ulsed  electron  Seam  experiments. 


6 .i+yPe  of  report  a period  covered 

■ Final  Repart  , „ 


T PER  FORM  i'ngTo  RG.  REPORT  NUMBER 


7.  author^; 

f 


8 CONTRACT  or  grant  NUMBERS 


I John  A.  /ilacFarlane,  Lt,  USAF 


T PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

AIR  FORCE  WEAPONS  LABORATORY 
Kirtland  AFB,  NM  87117 


10.  PROGRAM  ELEMENT  PROJECT.  TASK 
AREA  4 WORK  UNIT  NUMBERS 

62704H/WDNA4501 /Subtask 
N99QAXAA121,  Work  Unit  10 


II.  CONTROLLING  OFFICE  NAME  AND  AODRESS 


.12.  report  oate 

Novenber  T976 


>«yF»ets  ' ’ 


14  MONITORING  AGENCY  NAME  ft  ADDRESS^/  different  from  Controlling  Office) 


13.  TTOMggrfrfr  ! 

50  / 

IT  SECURITY  CL  AS577 ol  PrTs~ report) 


Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 


UNCLASSIFIED 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thru  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 

S'  J 

( / J 

J ^ 

— — 

18.  SUPPLEMENTARY  NOTES 

This  research  was  sponsored  by  the  Defense  Nuclear  Agency 

under  Subtask 

N99QAXAA1 21,  Work  Unit  10,  "Diagnostics  Development." 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


Heat  Transfer 
Calorimeter  Block 
Heat  Loss 


Diagnostics 
Electron  Beams 
Computer  Model 


Cooling  Curves 

Finite  Difference  Equations 


k 


ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

A study  to  analytically  model  the  thermal  behavior  of  a graphite  calorimeter 
block  subjected  to  an  initially  non-uniform  temperature  profile  has  been  con- 
ducted. The  model  allows  prediction  of  the  amount  of  energy  lost  from  the 
block  by  the  time  a uniform  temperature  profile  is  achieved.  It  is  assumed 
that  the  only  means  of  energy  loss  to  the  surroundings  is  that  of  thermal 
radiation;  the  effects  of  convection  and  thermal  conduction  to  the  surroundings 
are  neglected.  In  addition,  it  is  assumed  that  the  calorimeter  is  exposed  to.. 
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Figure  1.  Simplified  Model  of  Calorimeter 


Figure  2.  Typical  Temperature  History  Given  by  Thermocouple  Attached 
to  Rear  Surface  of  Calorimeter 


Shown  in  Figure  2 is  a typical  temperature  history  as  given  by  a 
thermocouple  attached  to  the  rear  surface  of  a calorimeter  block.  Because 
the  beam  energy  is  initially  deposited  in  a nonuniform  manner,  some  time 
must  elapse  before  it  is  uniformly  distributed  throughout  the  block  by 
means  of  thermal  conduction.  As  shown  in  Figure  2,  the  back  surface 
temperature  rapidly  rises  with  time  until  the  energy  is  uniformly  distributed, 
and  then  slowly  decreases  as  energy  is  lost  from  the  calorimeter  block  to 
the  surroundings. 

Because  the  thermocouple  measures  the  temperature  at  a given  point 
in  the  depth  of  the  block,  there  is  considerable  uncertainty  involved  in 
detecting  the  total  energy  deposited  unless  the  block  is  at  a uniform 
temperature.  For  the  case  presented  in  Figure  2,  it  is  generally  assumed 
that  the  block  is  at  a uniform  temperature  when  the  maximum  back  surface 
temperature  is  achieved.  The  total  energy  in  the  block  at  that  time  is 
then  obtained  using  the  specific  heat  of  graphite  at  that  temperature. 

To  obtain  the  energy  initially  deposited  in  the  calorimeter  block,  various 
extrapolations  of  the  "cooling  curve"  are  employed.  The  most  common 
extrapolation,  a linear  one,  is  shown  in  Figure  2.  This  extrapolation 
assumes  that  the  block  cools  at  a constant  rate  and  that  the  initial 
energy  in  the  block  can  be  given  by  the  initial  temperature  obtained  from 
such  an  extrapolation. 

It  is  the  purpose  of  this  study  to  analytically  model  the  thermal 
behavior  of  a graphite  calorimeter  block  subjected  to  an  initially  non- 
uniform  temperature  profile.  This  model  will  allow  prediction  of  the  amount 
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of  energy  lost  from  the  block  by  the  time  a uniform  temperature  profile 
is  achieved.  In  addition,  it  will  determine  the  validity  of  the  assumption 
that  the  block  is  at  a uniform  temperature  when  the  maximum  back  surface 
temperature  is  reached.  Then,  through  parametric  variations,  the  model 
will  be  used  to  determine  in  what  electron  beam  environments  graphite 
calorimetry  can  be  used  with  confidence. 

For  this  study,  the  simplified  calorimeter  model  of  Figure  1 is 
employed.  It  is  assumed  that  the  only  means  of  energy  loss  to  the  sur- 
roundings is  that  of  thermal  radiation;  the  effects  of  convection  and  thermal 
conduction  to  the  surroundings  are  neglected.  In  addition,  it  is  assumed 
that  the  calorimeter  is  exposed  to  a uniform  beam  with  an  area  equal  to 
that  of  the  calorimeter  face.  Thus,  conduction  in  the  lateral  directions 
is  neglected.  Finally,  it  is  assumed  that  the  beam  energy  is  deposited 
instantaneously,  with  the  resultant  temperature  profile  being  taken  as  an 
initial  condition. 

In  Section  II,  a mathematical  model  is  developed  for  a calorimeter 
block  subjected  to  these  conditions  and  constraints.  In  Section  III,  the 
model  is  applied  to  a graphite  block  through  specification  of  the  thermal 
properties  of  that  material.  Section  IV  describes  the  manner  in  which 
a computer  model  is  developed  to  determine  the  effectiveness  of  graphite 
calorimetry.  Section  V presents  the  results  of  a parametric  study  in  which 
the  radiation  losses  are  quantified  for  a variety  of  electron  beam  environ- 
ments. In  addition,  Section  V presents  an  evaluation  of  the  manner  in 
which  calorimetry  data  are  presently  reduced  to  yield  information  about 
the  loading  environment.  Finally,  Section  VI  presents  conclusions  and 
recommendations  arising  from  this  study. 
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SECTION  II 


DEVELOPMENT  OF  MATHEMATICAL  MODEL 


1.  INTRODUCTION 

For  the  problem  of  one-dimensional  heat  conduction  into  a graphite 
calorimeter  block,  proper  account  should  be  taken  of  the  variation  of 
thermal  conductivity  and  specific  heat  with  temperature.  For  this 
reason,  the  one-dimensional  heat  conduction  equation  cannot  be  reduced 
to  the  common  Laplace  equation.  Rather,  it  must  be  written  in  its  most 
general  form: 


_3T  _ 1_  3_  _3Ts 

3t  pC  3x  3x 


For  the  problem  being  addressed,  the  initial  and  boundary  conditions 
imposed  upon  equation  (1)  prevent  it  from  being  solved  through  analytical 
methods.  The  radiation  heat  transfer  at  the  free  surfaces  provides  a 
nonlinear  boundary  condition,  while  the  initial  temperature  profile  is 
generally  nonuniform.  For  these  reasons,  the  problem  being  considered  is 
best  analyzed  through  the  use  of  finite  difference  techniques. 

2.  DEVELOPMENT  OF  FINITE  DIFFERENCE  EQUATIONS 

Rather  than  apply  finite  difference  approximations  to  equation  (1), 
these  approximations  should  be  applied  to  the  First  Law  of  Thermodynamics 
for  the  calorimeter  block.  This  approach  allows  consideration  of  radiative 
losses  from  the  sides  of  the  block  without  resorting  to  a two-dimensional 
formulation. 
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Consider  the  calorimeter  block  shown  in  Figure  3.  If  this  block 
is  divided  into  M differential  elements,  each  Ax  thick,  an  energy  balance 
can  be  developed  for  each  element.  By  choosing  Ax  to  be  small,  the  ith 


element  can  be  considered  to  be  at  a uniform  temperature  T . Because  of 
geometrical  similarity,  the  development  will  be  identical  for  each  internal 
element.  To  account  for  the  front  and  rear  free  surfaces,  however,  a 
separate  development  is  necessary  for  the  1st  and  Mth  elements. 

a.  Internal  Elements 

Figure  4 illustrates  the  various  mechanisms  by  which  energy  is 
transferred  into  and  out  of  the  ith  internal  element.  Since  this  element 
is  at  a uniform  temperature  T , conduction  terms  for  the  y and  z directions 
can  be  neglected.  The  conduction  into  the  left  face  of  the  ith  element, 
given  by  results  from  a temperature  difference  between  element  i-1 

and  element  i.  Similarly,  Qj^x+£x)  represents  conduction  out  of  the  right 
face  of  the  ith  element.  Hie  energy  loss  to  the  surroundings,  assumed  to 

r 

occur  through  thermal  radiation  only,  is  given  by  Q^. 

The  First  Law  of  Thermodynamics,  as  written  for  the  ith  internal 
element , is 


3t 


k _ k r 

Qi(x)  " Qi(x+Ax)  ~ Qi 


(2) 


With  i denoting  position  and  n denoting  time,  the  left  hand  side  of 
equation  (2)  is  approximated  with  the  following  fully  implicit  finite 


difference  scheme: 


Note  that  the  subscript  n denotes  a known  quantity,  while  the  subscript 


n+1  denotes  an  unknown  quantity  at  the  next  time  step.  Using  the  same 


difference  scheme,  the  energy  conduction  terms  are  approximated  by 


i(x+Ax) 


3x  / (x+Ax) 


For  the  sake  of  simplicity,  without  any  loss  of  accuracy,  the  thermal  con 


ductivity  of  the  ith  element  is  calculated  using  the  known  temperature, 


The  radiation  loss  from  the  sides  of  the  element  is  given  by 


0 € 2Ax(Y+Z) 


Putting  equations  (3)  through  (6)  into  equation  (2)  and  rearranging 


yields  the  following  implicit  finite  difference  equation 
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Equation  (7)  is  a nonlinear  difference  equation  because  of  the  radiation 


term  involving  T.  . By  using  the  common  radiation  approximation  [1], 


4 3 

T,  ,,  = 4T.  T.  .. 

i,n+l  i,n  i,n+l 


3T. 


i,n 


(8) 


equation  (7)  can  be  linearized*. 


1.  Gaumer,  G.  R. , "Stability  of  Three  Finite  Difference  Methods  of 
Solving  Transient  Temperatures,"  ARS  Journal,  Vol.  32,  No.  10,  October  1962, 
pp  1595-1597. 

*The  radiation  approximation  is  obtained  by  taking  the  time  derivative 
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E. 
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Equations  (7)  and  (8)  can  be  combined  in  a considerably  simplified 
form  if  the  following  terms  are  defined: 


= 4 Ax  (Y+Z) 


YZ 


S„  . = 4ae  T.  ' 
2i  i,n 


3i 


S,  . 
4i 


. = ae["-3T.  4 - T 4] 

l |_  i,n  a J 


At 


2pC.  Ax 
i.n 


S,.  = k.  + k 
5i  i,n  i+l,n 


S£ . = k + 2k.  + k ... 

6i  l-l, n i,n  l+l, n 


S-,.  = k.  + k.  , 

7i  i,n  i-l,n 


(9) 


(10) 


(ID 


(12) 


(13) 


(14) 


(15) 


Each  of  these  quantities  are  known,  being  labeled  with  the  subscript  n. 
These  definitions  allow  equations  (7)  and  (8)  to  be  combined  in  the 
following  form: 
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T.  . (S,  . -S7.) 

l-l, n+1  4i  7i 


+ T.  ..  (-1  - S, . *S,.  - S. . -S  -S. .) 
i,n+l  4i  61  4i  1 2i 


+ Ti+1 , n+1  (S4i'S5i) 


= -T.  + (S.  . -S, -S,.) 

i,n  4i  1 3i 


(16) 


This  equation  is  the  final  form  of  the  finite  difference  approximation 
to  the  First  Law  of  Thermodynamics  for  the  ith  element.  It  can  be  used 
to  describe  each  of  the  internal  elements  depicted  in  Figure  3.  To 
completely  describe  the  calorimeter  block,  additional  equations  must  be 
developed  for  the  first  and  Mth  elements.  These  developments  are  presented 
in  the  following  subsection. 

b.  Elements  at  Front  and  Rear  Free  Surfaces 

Figure  5 illustrates  the  mechanisms  by  which  energy  is  transferred 
into  and  out  of  the  first  element.  This  element,  at  a uniform  temperature 
T^,  experiences  conduction  through  the  right  face  only.  Energy  transfer 
by  radiation  occurs  at  the  front  surface  as  well  as  from  the  free  side 
surfaces  of  the  element. 

The  First  Law  of  Thermodynamics,  as  written  for  the  first  element, 

I 

is 


• 

: 

I 

» 

: 
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1 k r 

3t  ^l(x=Ax)  ^1 


Using  the  same  finite  difference  approximations  employed  for  the  internal 
elements,  the  components  of  equation  (17)  are  given  by 


)E,  /3T,  \ f T,  - T. 

rr  - pVxYZW7  “ pCi,nAxYZ 


Ql(x=Ax)  - YZ 


(-k|i\  a ,YZ[v^2,„irh, 

V dx  y(x=Ax)  2 J 


n+1  Y2 , n+1 


= 0£  2Ax(Y+Z)  + YZ  T 
1 -L  > 


4 4 

- T 

n+1  a 


Putting  equations  (18)  through  (20)  into  equation  (17)  yields  a relation 
which  is  similar  to  equation  (7) 


T - T 

l,n+l  l,n 


- (k.  + k,  ) 

\ l>n  2,n  j 

+ ( k,  + k„  ) 

\ l,n  2,n^ 


2pC,  Ax 

l,n 


2 1 , n+1 


2pC,  Ax 
l.n 


2 2, n+1 


At  0 £ ( 4 Ax  (Y+Z) 

2pC . Ax2  V YZ 


i + 2YZ  Ax\ 

'■  / [ 1. 


4 4 

- T 
n+1  a 
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After  introducing  the  radiation  approximation  of  equation  (8),  the 
following  definitions  are  introduced  to  simplify  equation  (21) : 


S21  * 4°eTl,„' 


S3i  - ce 


— 3T  4 - T 4"| 
i l.n  a J 


2 

2pC,  Ax 
l.n 


Sri  = ki  „ + 

51  l,n  z,n 


Note  the  similarity  between  these  definitions  and  those  employed  for  the 
internal  elements.  Putting  equations  (22)  through  (25)  into  equation  (21) 
yields  the  following  final  form  of  the  finite  difference  equation  for 


the  first  element: 


- S4l'S51  - S4l'(Sl  + 24X)-S21> 


+T2,n+l(S41-S51> 


-I1,„+S41'<S1+24X)-S31 
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Using  a development  which  parallels  that  for  the  first  element,  the 
difference  equation  for  the  last  element  (i  = M)  can  be  written  as 


TM-l,n+l(S4M‘S7M^ 


+TM,n+l(  1 “ S4M-S7M  ~ S4M’(S1  + 2Ax)'S2M) 


-T„  + S,  • (S..  + 2 Ax)  • S„ 

M,n  4M  1 3M 


where 


S-,w  = 4oe.T„ 

2M  M,n 


S3M  = °£ 


j -3T  4 - T 41 

l M’n  a J 


54M  or  a 2 
2pCM,nAx 


S7M  = Sl.n  + ^Ti-l.n 


3.  SOLUTION  OF  FINITE  DIFFERENCE  EQUATIONS 

Equations  (16),  (26),  and  (27)  constitute  a set  of  M equations  which 
can  be  used  to  describe  the  thermal  behavior  of  the  calorimeter  block 
depicted  in  Figure  3.  Note  that  each  of  the  equations  is  of  the  form 
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Ki 


} a 

J v 


^ i 


for  i = M: 


AM=  s„.-s 


4M  7M 


(41) 


BM  ~1  S4M‘S7M  “ S4M*(S1  + 2Ax)‘S2M 


(42) 


Si 


= o 


(43) 


Dm  = -Tm  „ + S/M-(S1  + 2 Ax) • S 
M M,n  4M  1 3M 


(44) 


Because  equations  (16),  (26),  and  (27)  can  be  written  in  the  form  of 
equation  (32),  with  coefficients  described  by  equations  (33)  through  (44), 
they  constitute  a tridiagonal  matrix.  Use  of  these  equations  to  predict 
the  temperature  profile  history  for  the  calorimeter  block  requires  that 
this  matrix  be  solved  at  each  new  time  step.  The  matrix  is  easily  solved 
through  use  of  the  Thomas  Algorithm  [2] . 


For  each  new  time  step,  solution  for  T^  is  performed  as  follows: 


ei  = Bi  - 1 where  3i  = Bi 


(45) 


i-1 


D.  - A.y.  n 

i li-l  , 

Y±  = . where  y. 


6i 


(46) 


2.  Von  Rosenburg,  D.  U. , Methods  for  the  Numerical  Solution  of 
Partial  Differential  Equations,  pp  113,  American  Elsevier  Publishing 
Company,  Inc.,  New  York,  1969. 
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After  computing  6^  and  y^  for  all  i (called  the  First  Half  of  the  Thomas 
Algorithm),  one  then  computes  the  new  temperatures  by 


TM,n+l  YM 


r t 

T = y — * i+l>n+l 

i,n+l  Ti  g 


(called  the  Back  Half  of  the  Thomas  Algorithm) . 


4.  SUMMARY 

In  this  section,  the  First  Law  of  Thermodynamics  has  been  applied 
to  each  of  the  M difference  elements  represented  in  Figures  (3)  through 
(5).  This  has  resulted  in  a mathematical  model  which  is  to  be  used  to 
predict  the  thermal  behavior  of  graphite  calorimetry  under  a variety  of 
conditions.  Before  such  use  can  be  made  of  this  model,  information  is 
required  regarding  the  material  properties  of  graphite.  This  information 
is  given  in  the  following  section. 


SECTION  III 


MATERIAL  PROPERTIES  OF  GRAPHITE 


Before  the  mathematical  model  developed  in  Section  II  can  be  applied 
to  the  specific  case  of  a graphite  calorimeter,  values  must  be  assigned 
to  certain  material  properties.  These  properties  are  the  density,  p. 


the  emissivity,  e,  the  thermal  conductivity,  k,  and  the  specific  heat,  c. 

Graphite  is  produced  in  different  forms  through  various  means. 

For  this  reason,  different  values  of  density  are  quoted  in  the  literature. 

3 

The  value  used  in  this  study  is  1.8  g/cm  . 

The  total  hemispherical  emissivity  for  graphite  is  nominally  0.70 
to  0.95.  For  this  study,  a value  of  0.85  was  chosen. 

As  mentioned  in  Section  II,  the  thermal  conductivity  of  graphite 
is  a function  of  temperature.  For  the  range  of  temperatures  encountered 
by  graphite  calorimeters,  the  variation  of  the  thermal  conductivity  can 
be  quite  significant.  Reference  3 gives  values  of  k versus  T for  various 
forms  of  graphite;  for  each  case,  k was  measured  in  the  direction  of  the 
grain  axis.  A bilinear  relationship  was  found  to  well  describe  the 
behavoir  of  k,  which  varies  from  0.40  cal/ cm  sec  °K  at  0°K  to  0.05  cal/cm 
sec  °K  at  4000°K.  This  relationship  is  shown  in  Figure  6 and  is  represented 
mathematically  as 


0°K  < T < 1000°K  k = 0.400  - 2.50x10  T (cal/cm  sec  °K) 


3.  Goldsmith,  A.,  Waterman,  T.  E.,  Hirschnorn,  H.  J.,  Handbook  of 
Thermophysical  Properties  of  Solid  Materials,  Vol  I,  pp  115-121,  The 
MacMillan  Company,  New  York,  1961. 


Figure  6.  Specific  Heat  and  Thermal  Conductivity  of  Graphite  Versus  Temperature 
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1000°K  < T < 4000°K  k = 0.183  - 3.33x10  T (cal/cm  sec°K) 


The  specific  heat  of  graphite  is  also  a function  of  temperature, 
primarily  in  the  range  of  0°K  < T < 1000°K.  Reference  4 presents  a specific 
heat  versus  temperature  relationship  for  graphite.  Again,  a bilinear 
relationship  was  found  to  well  describe  the  thermal  behavior.  Shown 
in  Figure  6,  this  relationship  is 


0°K  < T < 1000°K  C = 5.0x10  T (cal/g°K) 


1000°K  < T < 4000°K  C = 0.50  (cal/g°K) 


Using  the  definition  of  specific  heat,  the  internal  energy  of  graphite 
can  be  expressed  as  a function  of  temperature  also: 


u(T)  = / C dT 

T 


where  u(T)  is  calculated  relative  to  some  arbitrary  temperature  T 


Taking  T to  be  the  ambient  temperature,  T , equations  (51)— (53)  yield 

O cl 

the  following  relationships  for  u(T): 


0°K  < T < 1000°K  u = 2.5xl0-4  (T2  - T 2)  (cal/g)  (54) 

— cl 


4.  Hultgren,  R.,  Selected  Values  of  the  Thermodynamic  Properties 
of  the  Elements,  American  Society  for  Metals,  Metals  Park,  OH,  p 91,  1973. 


21 


1000°K  <T<  4000°K  u = 0.50T  - 250(1.0  + 1.0x10  T ) (cal/g)  (55) 

" Si 

These  equations  can  be  written  in  terms  of  T versus  u: 

0 cal/g  1 u £ 228.5  cal/g  T = (4.0xl03  u + T 2)1/2  (°K)  (56) 

228.5  cal/g  < u T = 2u  + 500  (1  + lxlCf6T  2)  (°K)  (57) 

— a 

When  these  material  property  relationships  are  combined  with  the 
mathematical  model  of  Section  II,  the  result  is  a calculational  tool  for 
predicting  the  thermal  behavior  of  a graphite  calorimeter.  How  this 
calculational  tool  is  incorporated  into  a computer  model  for  performing 
such  predictions  is  the  subject  of  the  following  section. 


SECTION  IV 


COMPUTER  MODEL 

As  stated  earlier,  combining  the  mathematical  model  developed  in 
Section  II  with  the  material  property  relationships  defined  in  Section  III 
results  in  a calculational  tool  for  predicting  the  thermal  behavior  of  a 
graphite  calorimeter.  This  calculational  tool  has  been  incorporated  into 
a computer  model  to  allow  for  rapid  solution  of  a large  number  of  problems. 

The  computer  model  used  for  this  study  is  presented  in  the  following 
pages.  The  model  is  fairly  straightforward  and  should  be  self-explanatory. 
For  convenience,  all  symbols  employed  in  the  model  are  defined.  In  addition, 
various  sections  of  the  model  have  been  identified  with  regard  to  their 
function . 

The  following  section  presents  the  results  of  utilizing  this  model 
for  a variety  of  system  dimensions  and  initial  thermal  loading  conditions. 
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PROGRAM  MAC3(  IN»UT, OUTPUT> 


C 

C***  IOrNT IFTCATTON  OF  SYMBOLS  USED  IN  COMPUTE?  CODE 
C 

C •*  W - THICKNESS  OF  CALO=IMfTEx  9L0CK,  CM 
C **  Y - HEIGHT  OF  CAL0?IMET"F  BLOCK, CM 
C •»  1 - WID^H  OP  CALOPIMETPS  BLOCK, CM 

C **  PHO  - DENSITY  of  CALO=ImE'EP  block  vATEcIAL,G/CM**3 
C **  EMI  - TOTAL  ^EMISPHEPICAL  EMTSSINITY  OF  CALOFIMETEd  SURFACE 
C **  TA  - AM3IEN*  Tc MPE « ATU» E , K 

C **  SIG  - StcFAN  30L  T 7M  ANN  CONSTANT,  CAl/CM**?SEC  K**4 
C **  i - SUB$C=IpT  DENOTING  POSITION  IN  rALO=IMETEc  THICKNESS 
r **  x(I)  - LOCATION  of  ITh  ELEMENT,  CM 
C * * OX  - THICKNESS  of  each  DIFFERENCE  ELEMENT, cm 
o **  m - TOTAL  NUMBPP  of  DIFFrrENCE  ELEMENTS 
C **  N - SUBSCRIPT  OcNOTING  rIM? 

C **  TIM?  - ’’IMS  ELAPSED  SINCE  ENE°GY  DEPOSITION,  SEC 
C **  OT  - OUc  A ^ I ON  OP  ’"IME  STfP  ( \l  A P I A 3L  : ) , S " C 
C **  K - TOTAL  NUMBER  rF  tim'  STrPS  AlLOWEO 
C **  A3  - VARIABLE  USED  TO  ADJUST  NCXT  TtM£  step 
C **  A 3t  - (/  A 3 1 A PLF  USED  rr,  adjust  NEXT  TIME  STEP 
C **  Till)  - tFMPcpatUPE  OF  ITH  ELEMENT  AT  OLD  TIME  (N) 

C **  T 2 ( I ) - TPMO-PATUPP  QP  ITH  ELEMPNT  IT  NEW  TIME  <N*1) 

C.  **  UKI)  - SPPCIPIC  ENERGY  Op  ITh  flEM^NT  at  OLD  TIME  (N)  .CAL/GM 

C **  U2(I)  - SPECIFIC  ENERGY  OP  ITH  cLEMtNT  AT  NEW  TIME  (N*i>  , CAL/GM 

C **  C°(I)  - SPrCIPIC  HEA”  Or  ITH  ELEMENT, CAL/G-K 

C **  C<(I)  - THPJMAL  C ONCUCT I V I T Y DP  ITH  ELEMENT , CAL /CM* *2 -SEC-K 
C *»  UTOTl  - "OT„L  ENERGY  IN  CAL  OR I METE?  AT  OLD  TIME,CAL 
C **  UTDT2  - T QT  A L ENE=GY  IN  CALO’IMETER  AT  NCW  TIME, CAL 
C **  ELOSS  - PR R CrNT  EN-PGY  LOST  AT  cNO  OE  GIi/EN  TIME  STFP 
C **  SI  - FINITE  FIEP'PENCE  COEFFICIENT  DEFINED  IN  TEXT 

r **  S?(I)  TH ° DUG H S 7 ( I > - FINITE  DIFFERENCE  COPFFICIENTS  CEFINEC  IN  TEXT 
C **  A ( I ) tuc.ough  0(1)  - FI  N I T"  DIEFE?rNCP  C pf  Ff  IC I E NT  S DEFINED  IN  Tr  XT 
C **  BBCI).GG(I)  - v/APIABLES  USED  IN  SOLUTION  OE  M SIMULTANEOUS  EQUATIONS 
C 

DIMENSION  ”1(100) 

DIMENSION  T 2 ( 1 C 3 ) 

DIMENSION  U 1 ( 1 C 0 ) 

DIMENSION  U2 ( 1 0 3 ) 

DIMENSION  CP ( 1 C 3 ) 

DIMENSION  CKdCG) 

DIMENSION  X(100) 

DIMENSION  S2J1C0) 

DIMENSION  S ?( 1C  3 ) 

DI MrNS I ON  S4(100) 
dimension  :c(ir.o) 

DIMENSION  S€ ( 1 C 3 ) 

DIMENSION  S 7 ( 1 C 3 ) 

01  MR NS  I ON  A (10  0) 

DIMENSION  B(10C) 

OImpnSTCN  C (130 ) 

DIMENSION  C(10C) 

31 Mc N S I O N PR<13:) 

DIMENSION  GG  < 1 0 3 ) 

r. 

p***  d-Ftnjc  cyrr;,  0 I M E NS  I ON  S , '*  A T Ec  I A L r>POPEOTIES«  AND  APPROPRIATE  CONSTANTS 


L 


M7  = M/ 2 
9 X = W / M 
9T  = 1 .E-6 


C***S(7T  INITIAL  SPECIFIC  F.N"RGY  PROFILE 

data  Ul/  180  0..  1749. .1697  . , 164  6. , 1694.,  154 3. ,1491.,  144 C .,1389. ,1337 
...  1286.  . 1234.. 1133. ,1131. , 108L .. 10  29., 977. , 926. ,874. ,823. ,771,  , 
.72r..669.,617.,566.,514.,463.,411.,36C.,3C9..257.,206.,154., 

.1 J 3.  .51 . ,13.. 64*0 ./ 

C 

C***CALCULAT'r  INITIAL  Tf MPE°  ATUR E PPO^ILE  AND  TOTAL  ENERGY  FROM  SPECIFIC 
C*** FNER&Y  PROFILE 
UTOT  1=0  . 

09  4 1 = 1,  M 
X(I) = OX* ( I - 1 > 

1 IF (Ul (I > .GT  .228 . 5)G0  TO  2 

T1(I)=<4.0E3*U1(I)*TA**2)**.5 

GO  TO  3 

T1 (I) =2.*U1(I) *500 .*  < l.*l. 0E-6*T A**  2) 

UYOTl=UTQTi,ul (I)*RHO*OX*Y*7 
CONTINUE 


C***DPIMT  INITIAL  VALUES 

03 INT  F ,w, Y.7.PHO.EM1 ,TA.N,< 

5 R0P,1AT<1H1,/,1CX.*W=  * • t IQ  • 4,  * CM*,15X,*Y=  *.E10.4,*  CM*,14X. 

1*7=  *,<-10.4.*  C9*./,10x,*chO=  *.E  10. 4,*  G/CM3*.  10X,*EM1=  *.E10.4, 
215X,*Tr.  = * , E10 . 4 , * K*,/,1CX.*N=  * , 1 4 . 2 , * K=  *,I4,/7) 

PP I NT  6 

6 FOFNAT( 1CX, ’LISTING  OF  INITIAL  T-tfPcFAr  U E E AND  ENERGY  PROFILES*,/) 
PP  INT  7 

7 F9P'1AT(iax.2<lX,*I*,?X.*niSTANCE<CN>*,?X,*TF9PERATU°E<tO*,2X,*ENER 
*GYICAL/G)*.2X) ) 

09  9 1=1. M2 
I2=I*M? 

PRINT  8,I.X(I).Ti (I) « Ul ( I) ,I2,X(I2)«Tl(I2)  , U1 ( 12) 

FORMAT ( ( 10 x,2(I 3. 3X.F10 . 3. 3X.E11 .4, 6X, E1C . 4, 3X) ) ) 

CONT  INUE 


8 

9 

C 

C 


C*  * * O ALCUL  A T r VARIATION  OF  T -MPE  =>  A^UP  E 3 = OFILE  AS  A FUNCTION  Op  TIME 
C 

TIME=G . 

09  200  N = 1 , K 
TTMF=TIME*OT 
C 
C 
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C***OALCUL  ATE  THERMAL  CONO  UC  T I \l  I T Y AND  SPECIFIC  HEAT  FOP  EACH  NOCE 

no  12  1=1. m 

IF  <T1 ( I)  .GT  .1000  . > GO  TO  11 

CK(I)=.4C-2.5E-4*Ti(I) 

C.P(I)  =5.QE-4*T1(  I) 

GO  TO  12 

11  CK(I)=. 1833-3. 333E-5*T1(I) 

CP (I) =.30 

12  CONTINUE 
C 

C*  * * CALC  UL  A T r vAEIOUS  FINITC  OIFFEPENC"  COFFFICIENTS  OFFINEO  IN  TF  XT 
S1=(7+Y)*4.*(DX**?)/(Y*Z) 

S2(1)=4.*S:G*EM1*(T1(1)**3) 

S3(l)=(-3.*Tl(l)**4-TA**4)*3IG*EMi 
S4(1>=0T/(2.*PH0*CP(1) * ( OX  **2 ) » 

S3 (1) = OK (2) *CK ( 1) 

A ( 1 ) = 0 . 

9(1)=-1.-S4(1)*S9(1)-S4(1)*(S1*2.*DX)*S2(1) 

C( 1) = S4( 1) *S5( II 

0( U = -Tl (1) *S4 (1) *(2 .*OX«-Sl)*S3( 1) 

00  23  1=2.51 

S2C 11=4. *SI5*EM1* (T1 (I) **3> 

S3 ( I )= (- 3. *T1 ( I ) **4-TA**4) *SIG*EMl 
S*  ( I) = 0T/ ( 2.*PHO*CP( I) * (OX** 2) ) 

S5«H=CK(I  + l)*CK(I) 

S6(I)=CK(I+1)*2. *CK( I) *CK( I-l) 

S7  ( i ) =r.<  m ♦ck(I-i) 

A ( I) =S7 ( I) *S4( I) 

9(  I>=-1  .-S4(I)  *S9(I)  -S1*S4  (I>  *S2(I) 

C(I>=S4(U*S5(I) 
am=-Ti(u  ♦S4(i)*si*s3ii) 

20  CONTINUE 

S2  ( M) =4. *^10* EMI* (T1 (H) **?) 

S3 (M >=(-’. *T1 (M> **4-TA**4) *SIG*EMI 
54  < H ) =r<T/  (2  ,*RHO  *CP  C M ) * ( QX**2)  ) 

37(H) =OK (M) +CK (H-l) 
a(H)=S4(H)*S7(H) 

0(H)  = -i. -34(H)  *S7(H)  - S4('T)  *(S1*2.*0X)  *S2  (H) 

C(M) =0  . 

O(M) =-Tl (H) fS4(H) * (2.*0X*S  1)*S3(H) 

C 

c***ppoF0aM  SOLUTION  OF  M SIMULTANEOUS  EOUATIONS  9Y  THOMAS  ALGORITHM 

90(1) =n(i) 

GG  (i)  = r (i) /?(l) 

00  30  1=2. Ml 

19 (I)=3(I)-A(I)*C(I-i)7 99(1-1) 

GG(I)= (0(1) -A (I) *GG(I-1))/99(I) 

*0  CONTINUE 

99(H)  = =>(M)-A(H)*C(H- 1)  / 9a  ( M-l) 

GG(M)= (O(M) -A (H)*&G(H-1) ) /99(H) 

T 2 ( M)  - GG  ( M ) 

00  40  11=1, Ml 

1 = N - 1 1 

T2 (I) =GG (I) -(C ( I) *T2( 1*1) ) /39( I) 

40  CONTINUE 

C 

C***CALCUL ATr  SPECIFIC  PN:3&Y  PFQFILF  ANO  TOTAL  ENERGY  AT  NEW  TIME 


1 

I 


* 
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UT0T2=G. 
oo  5a  1 = 1, m 

IFCT’(I) .GT.1000 . ) GO  TO  4 2 
U2<I)=?.5R-4*(T2(I)**2-TA**2) 

GO  TO  4* 

42  'J2(I»  = .5*T?|I) -25  5.*  (l.*i,CE-6*TA**2) 

4 3 UTPT2  = U"r0T2+U2  (I)  *<?HO*OX*Y*Z 

50  CONTINUE 

ELOSS= (UTOT1-UTOT2)*10:./UTOT1 

c 

C***P®INT  OUTt>UT  FOR  SrLECTFO  TIMES 
9?  PRINT  93.TINr«N,0T 

98  FOFMAT(1H1,//.10X,*TIME=*.E1G,3,*  SEC*,5X,*N=  *,I5.5X,*DT=  *,E10.4 

..*  SEC*,//) 

P3 1 NT  99 

39  FORMAT ( 1CX, 2 (IX, * I * , 2 X , * 0 1 ST ANC E (CM) *, 2X, * 7 EMPER AT UP E ( <) * ,2X, *ENER 

*GY  <C AL/G) * . 2X)  ) 

00  111.  1=1, M2 
I3=I*M3 

p=>INT  iaO,I,X(I)  . T 2 ( I ) ,U2(I),I2,X(I2),T2(I2),U2(I2) 

IOC  FOFMAT ( <1CX,2(I3. 3X,E1C. 3, 3X.E1G ,4,6X,Elu  . 4,  3X)»  » 

1C1  CONTINUE 

PRINT  1C  2 , ELOSS 

102  FOkMAT(/,i;x.*PS3CFNT  ENERGY  LOST  = *,£10.4) 
r 

C*  * * ST  0D  PROGRAM  WHEN  T2(l)  REACHES  SPECIFIES  LE  rfEL 
105  T0IF=(T2 (1) -T A) 

IF(TOIF,LT.50.)GO  TO  330 
C 

C*  * * S I7F  OT  SAOJUSY  MAGNITUDE  OF  TIME  INCREMENT)  SUCH  THAT  MAXIMUM  TEMPERATURE 
C***CHANGP  WITHIN  9L0CK  IS  1 OEG°EP  K. 

A9  = C . 

00  13b  1 = 1.*' 

A3T=A°S(T2<I)”T1(I)) 

IF  I AOT , L T , A 3) GO  T0  1C6 
A3=  ART 

10b  CONTINUE 

OT=OT/ to 
C 

C***5PS=T  TEMPERATURE  ° ROF I L E FO°  NEW  TIME  STE3 

109  00  11C  1=1. M 

110  T11I)=T2(II 

2C;  CONTINUE 

c 

r 

300  CONTINUE 


t 


f- 

I 
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It  was  desired  to  investigate  the  thermal  behavior  of  graphite 
calorimetry  for  the  range  of  electron  beam  environments  encountered  in  material 
and  structural  response  testing.  To  accomplish  this,  the  calorimeter  system 
dimensions  were  fixed  and  the  developed  computer  model  was  exercised  for 
a variety  of  initial  thermal  loading  conditions.  The  calorimeter  block 
was  taken  to  be  square,  2 cm  x 2 cm,  and  1 cm  in  depth.  These  dimensions 
were  later  varied  over  a small  range  to  study  the  effect  of  such  variation. 

As  stated  in  the  Introduction,  the  initial  energy  deposition  profile 
is  dependent  upon  the  electron  beam  current  and  the  voltage  which  accelerates 
the  electrons.  Frequently,  however,  the  thermal  loading  conditions  are 
described  in  terms  of  mean  electron  beam  voltage  and  front  surface  specific 
energy  (cal/g).  For  this  study,  linear  energy  deposition  profiles  were 
chosen  to  represent  a range  of  electron  beam  environments.  A total  of 
twelve  profiles  were  used;  these  profiles  are  shown  in  Figure  7. 

Given  the  system  dimensions  and  an  initial  energy  deposition  profile, 
the  computer  model  provides  calculation  of  the  temperature  profile  history. 

An  example  of  such  a profile  history  is  shown  in  Figure  8.  The  calculations 
are  based  upon  a "1  MeV"  electron  beam  with  a front  surface  specific 
energy  of  1200  cal/g.  Through  thermal  conduction  into  the  graphite  block, 
the  energy  which  was  initially  deposited  in  a non-uniform  manner  quickly 
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becomes  more  uniformly  distributed.  During  this  equilibration  period, 
the  front  surface  temperature  is  rapidly  decreasing  while  that  of  the  rear 
surface  is  rapidly  increasing.  Soon  the  energy  becomes  uniformly  distributed, 
and  the  block  then  cools  slowly  through  radiative  losses  to  the  surroundings. 

For  the  example  of  Figure  8,  the  back  surface  temperature  reaches  a 
maximum  value  at  a time  of  3.07  seconds.  This  point  in  time  corresponds 
to  the  peak  in  the  cooling  curve  of  Figure  2.  As  discussed  in  Section  I, 
experimenters  generally  assume  that  the  calorimeter  block  is  at  a uniform 
temperature  when  this  maximum  back  surface  value  is  attained.  For  the 
case  presented  in  Figure  8,  this  is  seen  to  be  true.  This  trend  was 
observed  for  all  calculations  performed  under  this  study;  at  the  time  of 
maximum  back  surface  temperature,  the  maximum  and  minimum  temperatures  within 
the  block  always  differed  by  less  than  0.5  percent.  Thus,  for  the  range  of 
electron  beam  environments  presented  in  Figure  7,  the  assumption  of  a 
uniform  block  temperature  at  the  time  of  maximum  back  surface  temperature 
is  shown  to  be  valid.  This  is  considered  a significant  finding. 

Shown  in  Figure  9 are  continuous  temperature  histories  for  the  back 
surface  of  a calorimeter  block  subjected  to  a "500  keV"  beam.  For  front 
surface  specific  energies  of  600,  1200,  and  1800  cal/g,  the  peak  back 
surface  temperature  occurs  at  2.08,  2.19,  and  2.56  seconds,  respectively. 

Note  that  the  temperature  increases  at  the  same  rate  for  each  case,  but 

that  the  "cooling  curve"  becomes  steeper  as  the  front  surface  specific  ^ 

energy  is  increased.  Thus,  for  a given  mean  electron  energy,  one  would 

expect  that  as  the  front  surface  specific  energy  is  increased,  the  total 
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Figure  9.  Back  Surface  Temperature  Histories  for  a "500  keV"  Electron  Beam 


radiative  losses  at  the  time  of  maximum  back  surface  temperature  also 
increases . 

In  determining  the  total  energy  initially  deposited  in  the  calorimeter 
block,  various  extrapolations  of  the  cooling  curve  are  employed  by 
experimenters.  Given  a back  surface  thermocouple  measurement  as  shown  in 
Figure  9,  the  most  common  correction  involves  a linear  extrapolation  of 
the  cooling  curve  to  time  equal  zero.  The  temperature  obtained  by  this 
extrapolation  is  taken  to  be  that  which  would  result  if  the  calorimeter 
block  thermally  equilibrated  without  any  radiative  losses  to  the  surroundings. 
To  obtain  the  total  energy  deposited  in  the  block,  this  temperature  is 
multiplied  by  the  mass  of  the  block  and  the  specific  heat  of  graphite  at 
that  temperature.  The  linear  extrapolation  is  based  upon  the  assumption 
that  the  calorimeter  block  cools  at  a constant  rate.  Much  of  the  block, 
especially  the  front  surface,  actually  radiates  at  a decreasing  rate  since 
the  temperature  is  decreasing  with  time.  For  this  reason,  the  linear 
extrapolation  should  underestimate  the  amount  of  initially  deposited 
energy. 

To  quantify  the  error  involved  in  such  an  extrapolation,  the  amount 
of  energy  actually  lost  through  radiation  at  the  time  of  maximum  back  surface 
temperature  was  investigated.  These  results  are  shown  in  Figure  10  for 
each  of  the  environments  of  Figure  7.  It  is  seen  that  the  worst  case  for 
this  energy  loss  is  the  "1  MeV,"  1800  cal/g  electron  beam;  at  the  time  of 
peak  back  surface  temperature,  9.5  percent  of  the  initially  deposited 
energy  is  lost. 
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Shown  in  Figure  11  is  the  back  surface  temperature  history  for  the 
"1  MeV,"  1800  cal/g  electron  beam.  The  ■’.mount  of  energy  in  the  calorimeter 
block  at  the  time  of  peak  back  surface  temperature  is  given  by  that  tempera- 
ture to  be  522  cal.  Adjusting  this  value  for  the  amount  of  energy  lost 
through  radiation,  the  energy  initially  deposited  is  calculated  to  be 
577  cal.  If  a linear  extrapolation  of  the  cooling  curve  is  employed,  the 
temperature  obtained  yields  an  initially  deposited  energy  of  546  cal. 

This  amount  is  5.4  percent  below  the  correct  value.  Thus,  although  the 
linear  extrapolation  provides  a partial  correction  for  the  radiative  losses, 
it  does  indeed  underestimate  the  initial  energy  in  the  block.  Shown  in 
Figure  11  is  the  extrapolation  required  to  give  the  correct  value  for 
the  initial  energy. 

It  is  seen  that  for  the  cases  presented  in  Figure  10,  the  error 
involved  in  using  the  maximum  back  surface  temperature  to  calculate  the 
initially  deposited  energy  is  never  more  than  ten  percent.  Furthermore, 
a linear  extrapolation  of  the  cooling  curve  reduces  this  error  to  a maximum 
of  around  five  percent.  Based  upon  these  results,  it  is  concluded  that  for 
a calorimeter  of  dimensions  2 cm  x 2 cm  x 1 cm,  energy  lost  through 
radiation  is  not  a significant  problem,  and  a linear  extrapolation  of  the 
cooling  curve  is  an  adequate  means  of  accounting  for  such  energy  losses. 

To  investigate  the  effect  of  calorimeter  dimensions  on  the  amount  of 

2 2 

energy  lost,  the  calorimeter  area  was  varied  from  1 cm  to  25  cm  . It 
should  be  noted  that  to  maintain  conditions  for  one-dimensional  heat 
conduction  the  size  of  the  beam  was  always  taken  to  be  equal  to  the  area 
of  the  calorimeter.  Shown  in  Figure  12  are  the  results  of  this  parametric 
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Back  Surface  Temperature  (°K) 


Percent  Energy  Lost  (%) 


Figure  12.  Effect  of  Calorimeter  Area  Upon  Percent  Energy  Lost  at  Time  of 
Maximum  Back  Surface  Temperature 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

A computer  model  has  been  developed  which  provides  calculation  of  the 
temperature  profile  and  total  energy  within  a graphite  calorimeter  block 
subjected  to  rapid  energy  deposition  by  an  electron  beam.  The  model  is 
based  upon  an  assumption  of  one-dimensional  heat  conduction  with  radiative 
energy  transfer  occurring  at  all  free  surfaces  of  the  calorimeter  block. 

Through  a series  of  parametric  variations,  the  effectiveness  of 
conventional  graphite  calorimetry  has  been  investigated.  Results  indicate 
that  for  the  range  of  electron  beam  environments  typically  employed  in 
materials  and  structural  response  experiments,  the  problem  of  thermal 
radiation  is  not  significant.  For  the  worst  case  investigated,  only  ten 
percent  of  the  initially  deposited  energy  is  lost  before  the  calorimeter 
block  reaches  a uniform  temperature.  By  means  of  a linear  extrapolation 
of  the  cooling  curve,  the  initially  deposited  energy  can  be  computed  to 
within  five  percent  of  the  actual  value. 

The  developed  computer  model  can  be  used  to  further  reduce  the  error 
involved  in  estimating  the  total  deposited  energy.  Given  information 
regarding  the  diode  voltage  and  current  histories,  the  peak  back  surface 
temperature  measurement  can  be  utilized  to  yield  an  initial  energy 
deposition  profile  which  does  not  provide  any  correction  for  radiative 
energy  losses.  Using  this  deposition  profile  in  the  developed  computer 
code,  the  amount  of  energy  lost  through  radiation  at  the  time  of  maximum 
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back  surface  temperature  can  be  calculated.  Because  the  percent  energy 
lost  for  the  actual  deposition  profile  and  that  obtained  by  the  means 
described  above  will  not  differ  considerably,  the  deposited  energy  given 
by  the  peak  back  surface  temperature  need  only  be  adjusted  by  that  amount. 

These  conclusions  may  contradict  the  actual  experiences  of  some 
experimenters  . It  should  be  noted  that  the  developed  model  assumes  that 
the  only  mechanism  for  energy  loss  is  through  radiation  to  the  surroundings. 
In  real  experimental  configurations,  there  are  always  some  paths  for 
energy  to  leave  the  calorimeter  through  thermal  conduction.  An  effort 
should  be  made  towards  more  adequately  thermally  isolating  the  calorimeter 
system  from  its  surroundings.  If  this  could  be  done,  then  radiation  would 
be  the  only  energy  loss  mechanism.  This  study  has  shown  that  for  such 
conditions,  the  problem  of  thermal  radiation  can  be  adequately  accounted 
for. 
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LIST  OF  SYMBOLS 


C Specific  heat,  cal/g°K 

E Energy,  cal 

k Thermal  conductivity,  cal/cm  sec°K 

M Total  number  of  differential  elements 

L Time,  sec 

T Temperature,  °K  ^ 

Q Energy  flux,  cal/cm  sec 

W Thickness  of  calorimeter  block,  cm 

x Location  in  calorimeter  thickness,  cm 

Y Height  of  calorimeter  block,  cm 

Z Width  of  calorimeter  block,  cm 

Ax  Thickness  of  differential  element,  cm 

At  Duration  of  calculational  time  step,  sec 

e Total  hemispherical  emissivity  ^ 4 

0 Stefan  Boltzmann  constant,  cal/cm  sec°K 

p Density,  g/cm 

Factors  for  solution  of  simultaneous  equations 
Subscripts 

1 Evaluated  for  i th  differential  element 

n Evaluated  at  n th  time  step 

x At  the  location  x=x 

x+Ax  At  the  location  x=x+^x 
a Evaluated  at  ambient  conditions 

o Evaluated  at  some  arbitrary  value 

Superscripts 

k Conduction 

r Radiation 
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